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BLOCKED REGULAR FRACTIONAL FACTORIAL DESIGNS 
WITH MINIMUM ABERRATION^ 

By Hongquan Xu 

University of California, Los Angeles 

This paper considers the construction of minimum aberration 
(MA) blocked factorial designs. Based on coding theory, the concept 
of minimum moment aberration due to Xu [Statist. Sinica 13 (2003) 
691-708] for unblocked designs is extended to blocked designs. The 
coding theory approach studies designs in a row-wise fashion and 
therefore links blocked designs with nonregular and supersaturated 
designs. A lower bound on blocked wordlength pattern is established. 
It is shown that a blocked design has MA if it originates from an 
unblocked MA design and achieves the lower bound. It is also shown 
that a regular design can be partitioned into maximal blocks if and 
only if it contains a row without zeros. Sufficient conditions are given 
for constructing MA blocked designs from unblocked MA designs. 
The theory is then applied to construct MA blocked designs for all 
32 runs, 64 runs up to 32 factors, and all 81 runs with respect to four 
combined wordlength patterns. 

1. Introduction. Fractional factorial designs are widely used in scientific 
and industrial experiments. Blocking is an effective method for reducing 
systematic variations and therefore increasing precision of effect estimation. 
Experimenters often face the practical problem of choosing good fractional 
factorial designs and blocking schemes. 

Fractional factorial designs are typically chosen according to the minimum 
aberration (MA) criterion [12], which includes the maximum resolution cri- 
terion [1] as a special case. The study of blocking in fractional factorial 
designs is complicated by the presence of two defining contrast subgroups, 
one for defining the fraction and another for defining the blocking scheme, 
therefore, resulting in two types of wordlength patterns, one for treatment 
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and another for block. The MA criterion can be applied to the treatment 
and block wordlength patterns separately. However, MA designs with re- 
spect to one wordlength pattern may not have MA with respect to the other 
wordlength pattern. One approach, that taken by Sun, Wu and Chen [19] 
and Mukerjee and Wu [15], is to consider the concept of admissible blocking 
schemes, but it often leads to too many admissible designs. Another ap- 
proach is to combine the treatment and block wordlength patterns into one 
single wordlength pattern so that the criterion of MA can be applied to it in 
the usual way. Sitter, Chen and Feder [17], Chen and Cheng [5] and Cheng 
and Wu [10] have proposed four combined sequences, resulting in four MA 
criteria (to be defined later). See [9] for a related approach. 

A practical and important issue is how to construct MA blocked designs 
with respect to one or more criteria. This question is not adequately ad- 
dressed in the literature. Most of the existing MA blocked designs rely on 
the work of Sun, Wu and Chen [19], who obtained the complete catalog of 
blocked designs with 8, 16, 32, 64 and 128 runs for up to nine factors. MA cri- 
teria rank blocked designs according to the treatment and block wordlength 
patterns, which are often obtained by counting words in the treatment defin- 
ing contrast subgroups and alias sets. When the number of factors is large, 
there are a huge number of words to be counted, causing considerable dif- 
ficulties in computation. For example, when a design with 64 runs and 25 
factors is arranged in 8 blocks, there are 2^^ — 1 = 4,194,303 words to be 
counted. It is cumbersome and sometimes even impossible to do so for thou- 
sands or millions of different designs. This calls for alternative computational 
methods. 

To avoid the aforementioned computational difficulties, we take a coding 
theory approach and propose new methods to compare and rank blocked 
designs without using defining contrast subgroups and alias sets. The idea 
is originally due to Xu [22], who proposed the concept of minimum moment 
aberration and established its equivalence to MA for unblocked designs. We 
extend the concept of minimum moment aberration to blocked designs for 
three of the four MA criteria in Section 2. 

To further ease the computation burden, we study relationships among 
MA blocked designs under different criteria and develop a general theory 
on MA blocked designs. The coding theory approach studies designs in a 
row-wise fashion and therefore links blocked designs with nonregular and 
supersaturated designs. Results on nonregular and supersaturated designs 
are used to establish an important lower bound on blocked wordlength pat- 
tern. It is shown that a blocked design has MA with respect to all four 
criteria if it originates from an unblocked MA design and achieves the lower 
bound. It is also shown that a regular design can be partitioned into max- 
imal blocks if and only if it contains a row (i.e., treatment combination) 
without zeros. Sufficient conditions are given for constructing MA blocked 
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designs from unblocked MA designs. Some technical lemmas are presented 
in Section 3 and the main results are given in Section 4. We shall point out 
that, for simplicity, we focus entirely on regular designs, even though most 
of the results can be easily extended to nonregular designs. 

With the concept of minimum moment aberration and developed theory, 
we present methods to construct MA blocked designs in Section 5. We obtain 
MA blocked designs for all 8, 16, 27 and 32 runs, 64 runs up to 32 factors, 
and all 81 runs with respect to four combined wordlength patterns. The 
difference among MA blocked designs under different criteria is summarized. 

The rest of this section introduces some background. A regular s'^~^ de- 
sign is defined by k treatment defining words, which form the treatment 
defining contrast subgroup. The resolution [1] is the length of the shortest 
word in the treatment defining contrast subgroup. For i = 1, . . . ,n, let Ai^Q 
denote the number of words of length i in its treatment defining contrast 
subgroup. For two unblocked regular s"'~^ designs Di and D2, let r be the 
smallest integer such that ^^^0(^1) 7^ ^r,o(-C2)- Then Di is said to have less 
aberration than D2 if Arfi{Di) < Arfi{D2). If there is no design with less 
aberration than Di, then Di has MA. In short, the MA criterion sequentially 
minimizes ^1,0, ^2,0, • • • , ^n,o- 

To arrange a regular s"'~^ design in blocks of size s"'~^~p, one can 
choose p independent block defining words, which form the block defining 
contrast subgroup. There are (s^ — l)/(s — 1) block effects, each confounded 
with s'^ treatment effects. For i = 1, . . . ,n, let Ai^i denote the number of 
treatment words of length i that are confounded with some block effects. 

As done in the literature, we shall only consider regular main effect (RME) 
designs where none of the main effects is aliased with another main effect or 
confounded with a block effect. It is evident that, for RME designs, ^41^0 = 
^2,0 = ^1,1 = 0. The vectors Wt = (^3,0, • • • ,^n,o) and Wb = (^2,1, ■ ■ • , 
are called the treatment and block wordlength pattern, respectively. Let 
^0,1 = for convenience. 

MA criteria for blocked designs differ in how the treatment and block 
wordlength patterns are combined. Sitter, Chen and Feder [17] first proposed 
the combined wordlength pattern 

(1) iyscf = (^3,0,^2,1, Ao, ^3,1, A,0, ^4,1,- ••)> 

where Ai^i is ranked after for i = 2,3, Chen and Cheng [5] pointed 

out that the ordering of wordlength patterns in (1) violates the hierarchical 
assumption, and proposed the sequence 

(2) Wcc = (3^3,0 + ^2,1, Ao, 10^5,0 + ^3,1, ^6,0, • • ■), 

where the sum of (^*~^)^2«-i,o and Ai^i is ranked before ^2i,o for i = 2, 3, 

Cheng and Wu [10] proposed the two combined wordlength patterns 

(3) Wi = (^3,0, ^4,0, ^2,1, ^5,0, ^6,0, ^3,1, • • 
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(4) VF'2 = (^3,0, ^2,1,^0, ^5,0, ^3,l> ^6,0,- ••), 

where Ai^i is ranked after ^2j,o in Wi and after ^2i-i,o in W2 for i = 2, 3, 

We shall mention that sequence (4) was first proposed by Chen and Cheng [5] 
and later independently by Zhang and Park [25] and Cheng and Wu [10]. 

Four MA criteria result from sequentially minimizing the corresponding 
combined wordlength patterns. MA blocked designs under the W sequence 
are called MA W designs. 

An orthogonal array (OA) of A'" runs, n columns, s levels and strength t, 
denoted by OA{N,n, s,t), is an N x n matrix in which all possible s* level 
combinations appear equally often as rows for any set of t columns. 

2. A coding theory approach: minimum moment aberration. For a prime 
power s, let GF{s) be the finite field of s elements. Let Vn be the n- 
dimensional row vector space over GF{s), that is, Vn = {{vi, . . . ,Vn) -Vi € 
GF{s) for i = l,...,n}. 

An [n,m] linear code over GF{s) is a vector subspace of V„ with di- 
mension m so that it has distinct vectors. An [n,m] linear code D can 
be specified by an m x n generator matrix G whose rows form a basis for 
the code. Then D = {u GVn'-u = vG, v G Vm.}- A regular s"~'^ design is an 
[n,n — k] linear code over GF{s). For an introduction to coding theory, 
see [13], Chapter 4, and [20]. 

Consider arranging a regular s"'~'' design in equal-sized blocks. A design 
of this kind is called a regular (s"'"'^ : s^) design. Such a design is specified by 
a pair of matrices T and B, defined over GF{s) and of orders {n — k)xn and 
{n — k)x p, respectively, such that T has full row rank and B has full column 
rank. Then a typical block of the design consists of all level combinations 
of the form uT, with u € Vn-k and uB = v, where v is any fixed vector in 
Vp. Different blocks correspond to different choices of v. Since B has full 
column rank p, there are choices of v, leading to a division of the s"'~^ 
level combinations into blocks. See [9] and [15]. 

Let Lp = (s^ — 1) /{s — 1) throughout this paper. Suppose that the columns 
of B are bi,...,bp. Let F be the (n — k) x Lp matrix whose columns are 
A161 + • • • + Xpbp, where Aj G GF{s), at least one Aj / and the first nonzero 
Aj is 1. 

The columns of T and F can be viewed as points of PG{n — k — l,s), the 
projective geometry of dimension n — k — 1 over GF{s). In the terminology 
of projective geometry, F is a (p— l)-flat in PG{n — k — l,s). Then a regular 
(s""'' : s^) design is an RME design if and only if T and F are disjoint; see 
[5] and [15]. 

Let G = (T, F) be the (n — A;) x (n + Lp) matrix and D be the linear 
code generated by G. For convenience, write D = {Dt,Df), where Dt is 
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the N X n treatment matrix and Dp is the N x Lp block matrix, with 
N = s"'~^. For integers t > 0, define moments 

N N 

(5) Kt,o{D)=N-'J2J2MDT)]\ 

i=ij=i 

N N 

(6) Kt,,{D)=N"'Y.Y.MDT)r5,,{DF), 

i=ij=i 

where 5ij{DT) and 5ij{Dp) are the number of coincidences between the ith 
and jth rows of Dt and Dp, respectively. For two vectors u = {ui, . . . ,Un) 
and V = {vi, . . . ,Vn), the number of coincidences is the number of z's such 
that Ui = Vi. We take O'^ = 1 throughout the paper. 

Remark 1. The definitions of Ktfl{D) and Kt,i{D) given in (5) and (6) 
work for both regular and nonregular designs. For regular designs, the double 
summation can be replaced with a single summation; for example, (6) can 
be simplified to 

N 

(7) Kt,,{D) = N-'Y.i^i,{DT)]%,{DF), 

i=l 

where j can be any row number. 

Remark 2. Note that Dp is a rephcated OA{sP, Lp, s, 2). It follows from 
Lemma 1 of [14] that 6ij{Dp) takes on only two different values. Specially, 
let yi, . . . ,yN be the rows of Dp. Then 



(8) 5,,{Dt 



Lp = isP-l)/{s-l), iiyi = yj, 
Lp_i = Lp — s^"^, otherwise. 



For an integer k, let (^) = x{x - 1) ■ ■ ■ {x - k + 1) / k\ if k > 0, Q = 1 and 
(^) = if A; < 0. For integers k,j > 0, let S{k,j) be a Stirling number of the 
second kind, that is, the number of ways of partitioning a set of k elements 
into j nonempty sets. It is well known that S{k,j) = (1/j!) J2i=o{~^y~^ if)''''' 
for k> j > 0. For integers k,i>0, define 



k 

(9) Qk{i;n,s) = (-ir ^ j!5(fe,j>-^(s - 1)^-* 

For integers t,i>0, define 

k=0 



n — I 



(10) ct{i;n,s) = {s- 1) ^^(-l)'^ I n'~''Q,{z;n, s) 
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It is easy to show that S{k, k) = 1, Qk{k; n, s) = (— and Qk{i', n, s) = 
when i> k. Therefore, ct(t; n, s) = s~*(s — l)t! and ct{i; n,s) =0 when i>t. 

The following two lemmas regarding unblocked designs are from Xu [22, 
23]. 

Lemma 1. For a regular s'^~^ design D and integers t>0, 

min(t,n) 

(11) Kt,o{D)= ct{i;n,s)A^o{D), 

1=0 

where Ct{i;n,s) are constants defined in (10) and ylo,o(^) = l/(s — 1)- 

Lemma 2. Sequentially minimizing Ki^Q,K2fi, - ■ ■ ,Knfi is equivalent to 
sequentially minimizing Ai^O) ^2,0) • • • > ^n,o- 

The minimum moment aberration criterion [22] sequentially minimizes 
-^1,0) -^2, 0, • • • ) Knfl- Lemma 2 implies that the minimum moment aberration 
criterion is equivalent to the MA criterion for unblocked designs. 

Extending Lemma 1 to blocked designs, we have the following result. 

Theorem 1. For a regular : s^) design D and integers t>0, 

min{t,n) 

(12) Kt,i{D)=s'' ct{i;n,s)[A^i{D) + LpA^oiD)], 

i=0 

where Lp = (s^ — l)/{s — 1), ct{i;n,s) are constants defined in (10) and 
Ao,o{D) = l/{s-l). 

The proof of Theorem 1 requires the generalized Pless power moment 
identities, a fundamental result in coding theory. For clarity, all proofs are 
given in the Appendix. 

For an RME design D, Aifl{D) = A2fl{D) = Ao^i{D) = Ai^i{D) = 0. From 
(11) and (12), we obtain Kifi{D) = s'hi, K2fi{D) = s-^n{n + s-l), Kq^i{D) = 
s~^Lp and Ki^i{D) = s~'^nLp. Furthermore, 

(13) K\oiD) = 6s"^(s - l)A3fi{D) + s"^i(n2 + 3ns + s'^ - 3n - 3s + 2), 

(14) K2,i{D) = 2s"3(s _ l)A2,i{D) + s~Mn + s - l)Lp. 

We can define three minimum moment aberration criteria for blocked 
designs by replacing Ai^ and Ai^i with Kt^ and Ki^i in (1), (3) and (4). 
Because ct{t;n,s) is a positive constant, it follows from (11) and (12) that 
the minimum moment aberration criterion with respect to VKjcfj or W2 
is equivalent to its corresponding MA criterion. 
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The MA Wcc criterion defined in (2) is more complicated than the other 
three criteria. Nevertheless, from (13) and (14), we obtain 

Ks,o{D) + K2,i{D) = 2s-\s - l)[3A3,o(I?) + A2,i{D)] 

+ s"^n[(n^ + 3ns + - 3n - 3s + 2) + (n + s - l)Lp]. 

Therefore, minimizing K-^^ + -f^2,i is equivalent to minimizing 3yl3^o + ^2,i- 

3. Some lemmas. Suppose that D = {Dx,Df) is a regular (s"~'^ : s^) 
design. Let xi, . . . ,xn be the rows of Dt and yi,. . . ,yN be the rows of Dp, 
where N = s"'~^. For m = 1, . . . ,3^, let Dm be the s^-'^-p x n treatment 
matrix corresponding to the mth block. For integers t > 0, define moments 

gji — k — p gji — k — p 

i=i j=i 

where 5ij{Dm) is the number of coincidences between the ith and jth rows 
of Dm- Let Bm = {i : Xj is a row of Dm, 1 < i < N}. It is evident that i G Bm 
and j G Bm for some m if and only if yt = yj. It is useful to express Kt{Dm) 
in terms of the original design Dt as 

(15) K,{Dm) = s-'^^-'-P^ J2i^'^j(DT)t 

m 

Without loss of generality, assume that Di contains the null treatment 
(i.e., a row of zeros) and call Di the principal block. Then Di is an [n,n — 
k — p] linear code over GF{s) and other blocks Dm, 2 <m< s^, are cosets 
of Di; therefore, 

(16) Kt{Dm) = Kt{Di) for m = 2,...,sf. 

Note that D\ is possibly a supersaturated design in which the number of 
columns is larger than the number of rows. 

The next result shows that Kt^i{D) is determined by Ktfi{D) and Kt{Di). 

Lemma 3. Suppose that D is a regular {s^~^:s^) design and Di is its 
principal block. For integers t>0, Kt^i{D) = Lp^iKtfl{D) + s~^ Kt{Di) . 

Lemma 4. Suppose that D is an [s"'~^ : s^) RME design and Di is its 
principal block. Let J = n{s^~^~^~^ — 1){s"'~^~p — and rj be the frac- 
tional part of J. 

(i) Ki{Di) = s-'^n and Di is an OA{s"'-''-P,n,s,l). 

(ii) K2{Di) > s-{«-fc-p)[n2 + (s«-fc-P _ i)(j2 + ^(i _ j^))]. The equality 
holds if and only if the difference among all 5ij{Di), i <j, does not exceed 
one. 
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(iii) A'2,i(£>) >Lp_is-2^(n + s-l) + s^("-^'~P+i)[n2 + (s"-^-P-l)(j2 + 
77(1-77))]. 

A regular (s"~'^ : s^) design D = {Dt, Dp) can be viewed as an unblocked 
regular design. For clarity, denote this unblocked design as 

Liun- For integers t > 0, define moments Kt{D^n) = N'^EiLi J2jLi[Sij{Dnn)Y , 
where N = s"'~^ and 5ij{D^^) = 6ij{DT) + 6ij{Dp) is the number of coinci- 
dences between the ith and jth rows of -Dun- The next result shows that 
Kt{Dnn) is related to Ktfi{D), Kt~i^i{D), AV2,i(^), A't-2,o(-D) and so on. 

Lemma 5. For a regular (s"~'^ : s^) design D and integers t>0, 
Kt{D^^) = Kt,o{D) + tKt.i,i{D) 

r=2 ^ ' 

- LpLp_i{Lp~^ - LpZ\)Kt-r,o{D)]. 

The following two lemmas are useful to know when MA blocked designs 
are the same under different criteria. 

Lemma 6. If D has MA with respect to both Wsd and Wi, then D has 
MA with respect to W2- 

Lemma 7. Suppose there exists some constant < a < 3 such that 
a^3^0+^2,i is minimized for D. If D has MA with respect to both H^cf 
and W2, then D has MA with respect to Wcc- 

Lemma 7 is very useful to show the MA Wcc optimality. The condition 
a < 3 is necessary; see Section 5 for counterexamples. 

4. Main results. Lemma 4(iii) and (14) together yield a lower bound of 
^2,1 as follows. 

Theorem 2. For an (s""''' : s^) RME design D, 
A2,i{D)>[2{s-l)r^ 

x{-n{n + s - 1) + s-("-fc-P-2) 

x[n2 + (s"-'=-^'-l)(j2 + r?(l-7y))]}, 
where J = n{s^~^~^~^ — 1){s^~^~p — and rj is the fractional part of J. 



MINIMUM ABERRATION BLOCKED DESIGNS 



9 



Theorem 2 plays an important role in the theoretical development and 
construction of MA blocked designs. The lower bound is tight for p = n — 
k — 1 and n — k — 2. Note that an RME design achieving the lower bound 
does not always have MA. When s = 2 and p < n — k — 2, the lower bound 
can be improved in some cases if the results in [4] and [2] are used. However, 
the improvement is usually negligible, noting that ^2,1 must be an integer 
for RME designs. 

Corollary 1. If D has MA with respect to W^d and W2, and D 
achieves the lower bound in Theorem 2, then D has MA with respect to 

The next result provides a sufficient condition when MA blocked designs 
are the same under four criteria. 

Theorem 3. // Dt has MA among all regular designs and D 

achieves the lower bound in Theorem 2, then D has MA with respect to 
W,,i, Wi, W2 andW,,. 

When the lower bound in Theorem 2 is achieved, the principal block Di 
has minimum moment aberration among all x n designs. Theorem 3 

can be generalized as follows. 

Corollary 2. If Dt has MA among all regular s^~^ designs and the 
principal block Di has minimum moment aberration among all s"'~^~p x n 
designs, then D has MA with respect to Wgcf, Wi, W2 and Wcc- 

The next result gives a simple necessary and sufficient condition when a 
regular design can be partitioned into maximal blocks as an RME design. 

Theorem 4. A regular s^~^ design containing the null treatment can 
be partitioned into maximal blocks as an RME design if and only if 

it contains a row without zeros. 

Mukerjee and Wu [15] previously studied the maximal blocking problem 
with a projective geometric approach. They managed to obtain a complete 
solution for and designs. Our approach appears to be more pleas- 
ant than theirs. Theorem 4 gives a simple answer to the question. 

When s = 2, a row without zeros is necessarily a row of all I's. Then a 
row and its fold-over forms a block. The unblocked design must be a fold- 
over design. A regular fold-over design is also called an even design [11], 
because it contains only words of even length. Whether or not a design is an 
even design can be simply checked by its wordlength pattern. The following 
corollary is a special case of Theorem 4. 
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Corollary 3. A regular 2""'^ design containing the null treatment can 
be partitioned into maximal 2"^"^"^ blocks as an RME design if and only if 
it is an even design. 

It is of special interest to know when an unblocked MA design can be 
partitioned into maximal blocks. Unblocked MA 2'^~^ designs were given 
by Chen and Wu [8] for A: = 1,2,3,4 and by Chen [6] for k = 5. Combining 
their results and Corollary 3, we have the following result. An unblocked 
MA 2"-^ desi gn can be partitioned into maximal 2" ^ blocks as an RME 
design as follows: 

(i) when k = 1 and n is even, 

(ii) when k = 2 and n is a multiple of 3, 

(iii) when k = 3 and n = 7t + q for integers t > and g = 7, 11, 

(iv) when A; = 4 and n = 15t + q for integers t > and q = 8, 12, 15, 20, 

(v) when k = 5 and n = 31t + q for integers t>0 and q = 16, 21, 24, 28, 31, 
37,40,44. 

Furthermore, it is known from coding theory that even designs are the 
only designs of resolution IV for 5A^/16 <n < N/2 with N = 2""^^; see [3]. 
For such n, an unblocked MA 2"'~^ design can always be partitioned into 
maximal 2"'~^~^ blocks as an RME design. 

To describe the next result, let F he an {n — k — 2)-flat and T be the 
complement of F in PG{n — k — l,s). Let H^-k be the linear code generated 
by T. Note that H^-k is unique up to isomorphism. It is evident that Hn-k 
and its projection designs (i.e., subsets of columns) can be partitioned into 
maximal s"'~'^~^ blocks as RME designs. The reverse is also true in the 
following sense. If an s"'~^ design can be partitioned into s"'~^~^ blocks as 
an RME design, then it is isomorphic to a projection design of Hn-k- The 
next result characterizes MA (s"-'^ ; §"-'^-1) RME designs. 

Theorem 5. If Dt has MA among all projection designs of Hn-k, then 
Dt can be partitioned as an {s"'~^ : s"""^"^) RME design D that has MA 
with respect to Wscf, Wi, W2 and Wcc- 

Theorem 5 shows that MA blocked (s""'^ : s""^'""*^) designs are the same 
for all four criteria when they exist. As a numeric illustration, consider 
s = 3. For 27 runs, H^ is the unique MA 3^~^ design. According to Xu 
[23], for 4 < n < 9, MA 3'^-("-3) dggig^g are projection designs of the MA 
3^~^ design; therefore, they can be partitioned into maximal 9 blocks and 
resulting RME designs have MA with respect to all four criteria. For 81 runs, 
Hi is the unique MA 327-23 (jggigji^ According to Xu [23], for 5 < n < 9 
and 12 < n < 27, MA 3"~("-4) designs are projection designs of the MA 
227-23 design; therefore, they can be partitioned into maximal 27 blocks 
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and resulting RME designs have MA with respect to all four criteria. For 
n = 10, 11, MA 3"-("-4) dggigj^g are not projection designs of the MA 327-23 
design; therefore, they cannot be partitioned as RME designs with 27 blocks. 
The second best designs are projection designs of the MA 327-23 figgig^; 
therefore, they can be partitioned into maximal 27 blocks and the resulting 
RME designs have MA with respect to all four criteria. 

When s = 2, Hn~k is an even design with resolution IV. We have the 
following result. 

Corollary 4. If a regular 2'^~^ design has MA among all even designs, 
then it can be partitioned as a 

^^n-k . 2n-fc-i) ^j^^^ design that has MA with 
respect to Wgd, Wi, W2 and Wcc- 

Recall that a regular (s"'"'^ : s^) design D can be viewed as an unblocked 
regular gi'i^+Lp)-ik+Lp) ^jggjgj^ Dun- The next result provides a sufficient con- 
dition when an MA blocked design originates from an unblocked MA design. 

Theorem 6. // Dx has MA among all regular s"'~^ designs and the 
unblocked design Dun has MA among all regular s(^+Lp)-{k+Lp) d^g-ig^s, then 
the blocked (s"""^ : s^) RME design D has MA with respect to Wsc{, Wi, W2 
and Wcc- 

Theorem 6 is most useful when p = l. It happens frequently that an un- 
blocked MA s""'^ design can be extended to an unblocked MA 
design by adding an extra column. For example, according to Chen, Sun 
and Wu [7] and Xu [23], for 8 and 27 runs, MA unblocked designs are in 
sequential order for all n. Whenever this happens, the extra column can be 
used as the block generator, and the resulting [s^~^ : s^) design has MA with 
respect to VFscfj W2 and Wcc- 

Theorem 6 is less useful when p> 1 because Dun usually does not have 
MA. The following result is interesting in this regard. 

Theorem 7. // Dt has MA among all regular s"~^ designs and D has 
MA with respect to Wgd, then D has MA with respect to both Wi and W2- 
If, in addition, aA3^Q{D) + A2^i{D) is also minimized for some constant 
< a < 3, then D has MA with respect to Wcc- 

Theorem 7 implies that when MA blocked designs are different under 
Wsc{, Wi and W2-, an MA W^d blocked design must not originate from an 
unblocked MA design. 
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5. MA blocked designs. MA blocked designs with respect to T^^cf, Wi 
or W2 can be obtained by computing and comparing moments Kt^Q and Kt^i 
for all possible blocking schemes. This is a feasible task when the number 
of blocking schemes is not too large; see [24] for details, where MA blocked 
designs for all 32 runs, 64 runs up to 32 factors, and all 81 runs with respect 
to Wsci, Wi and W2 are given. 

However, this method cannot be used to construct MA Wcc designs be- 
cause there is no equivalent minimum moment aberration criterion with 
respect to (2). Furthermore, an essential difference exists between the MA 
Wcc criterion and the other three criteria. Because the MA Wsd , Wi and W2 
criteria minimize ^3^0 first, there is no need to search over resolution III de- 
signs whenever blocking schemes from resolution IV designs exist. However, 
the MA Wcc criterion minimizes 3^3^0 + ^2,1 first. Combining ^3^0 with ^2,1 
makes it more difficult to construct MA Wcc designs than other types of MA 
designs. To determine the minimum of 3j43^o + ^2,1 1 a simple strategy is to 
search over all resolution III designs. This requires a complete catalog of 
resolution III designs, but such a catalog is not available for 64-run designs. 

Combining the developed theory and computer search, we obtain MA Wcc 
designs for all 8, 16, 27 and 32 runs, 64 runs up to 32 factors, and all 81 
runs. Previously, Chen and Cheng [5] developed a theory to characterize MA 
Wcc designs in terms of their blocked residual designs and obtained MA Wcc 
designs for all 8 and 16 runs and 32 runs up to 20 factors. 

Here we explain how to construct MA Wcc designs for 64 runs and n < 32 
with the results of Xu and Lau [24]. First, for p = 5 and 6 < n < 32, by 
Theorem 5, MA (2'^-("-6) : 2^) desi gns are the same under all four criteria; 
therefore, MA designs given by Xu and Lau [24] have MA with respect to all 
four criteria. Indeed, they can be easily constructed by searching over MA 
projection designs of the unique even 2^^~^^ design. Next, for p = 1, because 
unblocked MA 2"-("-'') desi gns are in sequential order for n = 6-7, 8-12, 
14-20 and 21-33, by Theorem 6, we obtain MA (2"-(""6) : 2^) designs with 
respect to all four criteria for all 6 < n < 32 but n = 7, 12, 13 and 20. For 
n = 7, 12, 13 or 20, according to [24], MA W2 and Wgcf designs coincide and 
have ^2,1 = 0; therefore, by Lemma 7, MA Wcc designs also coincide with 
MA W2 and Wscf designs. 

The situation for p = 2, 3, 4 is more complicated than that for p = l and 
5. We first compute the lower bounds of ^2,1 in Theorem 2, which are 
given in Table 1. It is evident that a lower bound can be replaced by the 
smallest nonnegative integer that exceeds it if it is negative or not an integer. 
According to Xu and Lau [24], MA W2 designs achieve the modified lower 
bounds of ^2,1 except for the following 22 cases: p = 2, n = 19-26, 31, 32; 
p = 3, n = 29-32; and p = 4, n = 25-32. Furthermore, MA Wscf and W2 
designs coincide except for p = 2 and n = 7, 12. Then, by Corollary 1, MA 
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Table 1 

Lower bound of A2,i in Theorem 2 for 64 runs, 6 < n < 32 and p = 2, 3, 4 



p 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


2 


-1.5 


-1.5 


-1.5 


-1.5 


-1.2 


-1.3 


-0.8 


-0.7 








1 


1.2 


2.3 


2.7 


3 








1 


1.5 


2.5 


3.5 


4.5 


6 


7 


9 


10.5 


12.5 


14.5 


16.5 


4 


3 


5 


7 


9 


12 


15 


18 


22 


26 


30 


35 


40 


45 


51 


P 


20 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 


31 


32 




2 


3.8 


4.5 


5.5 


6.5 


7.5 


8.7 


9.8 


11.3 


12.2 


14 


15 


17 


18.2 




3 


19 


21 


24 


26.5 


29.5 


32.5 


35.5 


39 


42 


46 


49.5 


53.5 


57.5 




4 


57 


63 


70 


77 


84 


92 


100 


108 


117 


126 


135 


145 


155 





Wsci and W2 designs also have MA with respect to W^c except for the 24 
special cases, which require additional computer search. 

For the 24 special cases, Theorem 2 and Lemma 7 are again used to ease 
computation. Consider, for example, p = 4 and n = 29. According to Xu 
and Lau [24], MA Vl^scf) W^i and W2 designs coincide and have ^3^0 = and 
^2,1 = 196. The lower bound of A2^i is 126. To determine the minimum of 
3^43^0 + ^2,1, we only need search over all designs with ^3^0 < (196 — 126)/3 = 
23.3, leading to A^.^ < 23. This is a feasible task. A complete enumeration 
(to be explained later) shows that there are exactly 17 regular 2^^~^^ designs 
with ^3^0 ^ 23, among which one has resolution IV. It is straightforward to 
verify that 2.9^3^0 + ^2,1 has minimum 196 among all 17 2^^~^^ designs with 
^3,0 ^ 23. Then, by Lemma 7, MA Wcc designs coincide with MA W2 and 
Wscf designs. 

When MA W2 and Wgcf designs are different or when they do not minimize 
aAs^o + ^2,1 for all a with < a < 3, Lemma 7 cannot be used; then MA Wcc 
designs are determined by sequentially comparing the complete sequence 
in (2). Fortunately, this happens only for the following five cases: {n,p) = 
(7, 2), (12, 2), (25, 4), (26, 4), (29,3). For the first two cases, MA Wcc designs 
coincide with MA W2 designs; for the last three cases, MA Wcc designs 
are different from MA W2 designs, which coincide with MA Wsd and Wi 
designs. 

Now we explain how to enumerate all 2^^~^^ designs with yl3^o ^ 23. Note 
that a 3-letter word consists of three factors and there are 29 factors in a 
229-23 (jggign. Therefore, for any 2^^~^^ design with ^3^0 = 23, there must 
exist a column appearing in at least 3 x 23/29 = 2.4 or 3 words of length 3. 
Deleting that column yields a 2^^~^^ design with A3 < 20. Therefore, all 
229-23 designs with ^3^0 ^ 23 can be enumerated by adding a column to all 
228-22 (jggigjjg with ^3^0 ^ 20, which in turn can be enumerated by adding a 
column to all 2^''"^^ designs with ^3^0 ^ 17- This can be done sequentially in 
the same way as in Chen, Sun and Wu [7] and Xu [23], as long as the number 
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of designs is not too large at each step. We shall point out the importance of 
the lower bound of ^2,1 in Theorem 2. Without this bound, one has to search 
over all 2^^~^^ designs with ^3^0 ^ 196/3 = 65.3. This is not a feasible task 
because there are more than 100,000 2^^~'^^ designs with ^3^0 ^ 65 and it is 
impossible to enumerate all of them with the current method and computer. 

Finally, we summarize the differences of MA blocked designs under differ- 
ent criteria for all 8, 16, 27, 32 runs, 64 runs up to 32 factors, and all 81 runs. 
We observed that MA blocked designs under all four criteria are the same 
in most cases. This occurs for all 8 and 27 runs, which can be easily verified 
with Theorems 5 and 6. When MA blocked designs under four criteria are 
not all the same, one of the following four situations occurs: 

1. MA Wi, W2 and Wcc designs are the same, but they differ from MA Wsd 
designs. 

2. MA W2, Wsc{ and Wcc designs are the same, but they differ from MA Wi 
designs. 

3. MA Wi, W2 and H^cf designs are the same, but they differ from MA Wcc 
designs. 

4. MA W2 and Wcc designs are the same, but they differ from MA Wi or 
Wsci designs. 

Situation 1 occurs once for 32 runs with {n,p) = (6, 1) and once for 64 runs 
with {n,p) = (7,2), and does not occur for 16 and 81 runs. Situation 2 
occurs twice for 16 runs with (n,p) = (5,1), (5,2), 12 times for 32 runs, 
35 times for 64 runs, and twice for 81 runs with {n,p) = (11, 1), (11, 2). 
Situation 3 occurs once for 32 runs with {n,p) = (13,3), three times for 
64 runs with {n,p) = (25,4), (26,4), (29,3), and three times for 81 runs with 
{n,p) = (9, 2), (17, 2), (21, 2). Situation 4 occurs only once for 64 runs with 
(n,p) = (12,2). 

Except for situation 3, MA Wcc designs coincide with MA W2 designs, 
which are given by Xu and Lau [24]. Table 2 gives MA designs for situa- 
tion 3 with treatment and block columns in the same fashion as Cheng and 
Wu [10] and Xu and Lau [24] . The designs are labeled as n-k • i /Bp{W) , where 
i denotes the rank of the unblocked s"'~^ design under the MA criterion, p 
denotes the number of block variables, and W denotes the MA ly-criterion. 
See Xu and Lau [24] for generator matrices and column labels. To save space, 
in Table 2 independent columns are omitted in the treatment columns; only 
generators are given in the block columns, treatment wordlength pattern 
is truncated as Wt = (^3^0; ^4,0) ^5,0; ^6,0) and block wordlength pattern is 
truncated as Wi, = (^2,1, ^3,1, ^4,1, ^5,1). The last two columns in Table 2 
give the numbers of clear main effects (CI) and of clear two- factor interac- 
tions (C2). A main effect or two-factor interaction is clear if it is not aliased 
with any other main effect or two-factor interaction and is not confounded 
with any block effect [19]. 
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Table 2 

MA blocked designs for situation 3 



Design Treatment Wt Block Wb C1C2 



32 runs 



13-8.1/B3(*) 


31 7 11 21 25 13 14 19 


55 96 


3 5 17 


36 310 


13 





13-8.4/B3(V[/cc) 


31 7 11 21 13 14 26 3 


4 39 32 48 


5 10 19 


22 76 124 288 


4 







64 runs 










25-19. 1/B4(*) 


31 35 13 52 14 55 37 61 11 19 435 


3 5 9 48 


144 5923 


25 







21 44 7 62 25 49 22 41 38 


5440 










25-19. 17/B4(W^cc 


) 31 35 13 52 14 55 21 37 11 19 8 378 336 


3 5 17 41 92 568 2688 


8 







25 38 7 26 49 22 28 50 9 


4032 




13104 






26-20. 1/B4(*) 


31 35 13 52 14 55 37 61 11 19 515 


3 5 9 48 


156 6999 


26 







21 44 7 62 25 49 22 41 38 26 


7062 










26-20.50/B4(M^cc 


) 31 35 13 52 14 55 21 37 11 19 16 386 672 


3 5 17 41 100 632 3248 










25 38 7 26 49 22 28 50 9 33 


4368 




15960 






29-23. 1/B3(*) 


31 35 13 52 14 55 37 61 11 19 819 


5 17 33 


91 5187 


29 







21 44 7 62 25 49 22 41 38 26 


14560 












28 42 47 












29-23.4/B3(M^c:c) 


31 35 13 52 14 55 37 61 11 19 12 707 640 


9 20 38 


46 484 2252 


4 







21 44 7 62 25 49 22 41 26 28 


11536 




14016 








42 56 3 














81 runs 










9-5.1/B2(*) 


22 9 24 31 34 


18 36 12 


4 20 


9 30 117 162 


9 





9-5.2/B2(W^cc) 


22 9 24 31 3 


1 18 27 28 


6 18 


6 44 90 186 


6 


5 


17-13. 1/B2(*) 


22 9 24 31 3 25 13 37 6 18 


20 336 1014 4 15 


40 210 2079 










7 35 12 


5072 




9256 






17-13.2/B2(Wcc) 


22 9 24 31 3 25 13 37 6 18 


23 306 110712 15 


28 303 1782 










7 35 16 


4952 




9814 






21-17.1/B2(*) 


22 9 24 31 3 25 13 37 6 18 


51 729 3717 4 26 


48 550 4590 










7 35 12 38 15 16 19 


21819 




32418 






21-17.2/B2(W^cc) 


22 9 24 31 3 25 13 37 15 23 


52 720 3735 11 30 


45 573 4545 










16 34 6 38 7 18 26 


21876 




32310 







Notes. (*) MA Wsd, Wi and W2 designs. 



For all designs given in Table 2, MA Wcc designs have a larger ^3^0 value 
but a smaller A2^i value than corresponding MA designs under the other 
three criteria. Indeed, these MA Wcc designs achieve the lower bound of ^2,1 
in Theorem 2, whereas MA designs under the other three criteria originate 
from unblocked MA designs. 

Note that for 81 runs, MA Wsc{,Wi and W2 designs 9-5.1/B2 have the 
same 8^3,0 + ^2,1 value as the MA Wcc design 9-5.2/B2. This also happens 
with 21-17.1/B2 and 21-17.2/B2. These examples show that the condition 
a < 3 in Lemma 7 is necessary. 
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APPENDIX 

Further notation and results in coding theory are necessary in order to 
prove Theorem 1. The Hamming weight of a vector u = (ui, . . . , denoted 
by wt{u), is the number of its nonzero components Uj. 

Associated with any hnear code D is another hnear code, called its dual 
and denoted by D-^. Suppose D is an [n,m] linear code with generator 
matrix G over GF{s). The dual D-^ is the null space of G, that is, D-^ = 
{u ^Vn'- uG' = 0}, where G' is the transpose of G. The dual is indeed 
the defining contrast subgroup of D. 

Suppose D is an [ni + n2,m] linear code over GF{s) and is its dual 
code. Each vector n in D and D-^ can be written as n = {ui,U2), where 
ui € and U2 G Vn2 - Let Bi^^i2{D) and Bi^^i2{D^) be the number of vectors 
in D and respectively in with wt{ui) = ii and wt{u2) = i2- 

The following result, a special case of Lemma 4.3 of Xu [21], generalizes 
the Pless power moment identities [16]. 

Lemma A.l. For integers ki,k2> 0, 

ni n2 ni 712 

^-m ^ g ilH'^'B,,,,,{D) = E E B,,,,,{D^)Q,,{ji;ni,s)Qk,{j2;n2,s), 

n=0i2=0 ji=0j2=0 

where (5fc(J;^^,s) is defined in (9). 

Proof of Theorem 1. Let N = s"-~'' and n2 = Lp. Then D = {Dt, Dp) 
is an [n + n2, n — k] linear code. Let be the dual code of D. Each vector u 
in D and D-^ can be written as u = {ui,U2), where ui G and U2 € Ki2- It 
is known that the wordlength patterns are proportional to the split weight 
distributions of as follows: for i = 0, . . . , n, 

(A.l) Aifl{D)=B,,o{D^)/{s-l) and = - 1); 

see [18] and [5]. By (7), 

n n2 

Kt,i{D)=N'' Y.in-ii)\n2-i2)B^,,,2{D) 

il=0j2=0 
n n2 t / , \ 

= E E E U i-l)''n'-'^'^{n2 - i2)i?n.2p) 
ii=0i2=0fc=0 ^ ^ 

= E u (-1)'^*^' E E - i^2)i?n,i2(^)- 

/c=0 ^ ^ 41=012=0 



MINIMUM ABERRATION BLOCKED DESIGNS 17 

By Lemma A.l, 

KtAD) = E U (-l)'n*-'^ E E B,,,,{D^)Q,{jr,n,s) 

X [Qo(i2;?^2,s)n2 - (5l(J2;^^2,s)]• 
Recall that Qkij',n, s) = for j > k. Then 

Kt,i{D) = (I) (-l)V-^ J2 Qk{ji;n,s)A{D^,n;n2,s), 
k=0 ^ ^ ji=0 

where 1^(1)^ , 2\\n2, s) = Bj^fi{D^)Qo{0;n2, s)n2 - Bj^^o{D-^)Qi{0;n2, s) - 
Bj-^^i{D^)Qi{l;n2, s). Note that <5o(0;n2,s) = 1, Qi(0;n2,s) = ^2(5 - l)s~^ 
and Qi(l; n2, s) = —s~^. Then 

fc=0 ^ ji=0 

X [B,,,o(^^)n2S-' + Bj,,^{D^)s-^] 
= j:i:(l){-l)'n'-'QkUi;n,s) 

ji=Ofc=0 ^ ^ 
n 

= E ctUi;n,s)is-l)-^[B,,,o{D^)n2 + B,,^i{D^)]s-\ 

Then (12) follows from (A.l) and the fact that ct{ji;n,s) = when ji > t. 
□ 

Proof of Lemma 3. By (8) and (15), 

N N 

i=ij=i 

+ ^''E E E mDT)]\L,-L,^,) 

sP 

= Lp^iKtfl{D) + (Lp - Lp_i)s-2P ^ Kt{D^). 

771=1 

Then the result follows from (16). □ 
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Proof of Lemma 4. Let Ni = s''-''~p and Jt{Di) = Y.i<i<j<N, [%(^i)]V 
[Ni{Ni - l)/2] for t > 0. It is easy to verify that, for t > 0, 

(A.2) Kt{Di) = N{'[{Ni - l)Jt{Di)+n% 

(i) Recall that for an (s""^ : s^) RME design D, Kifl{D) = s'^n and 
Ki i{D) = s-^nLp. Then, by Lemma 3, Ki{Di) = s[Ki i{D)-Lp^iKi o{D)] = 
s-^n. By (A.2), MD,) = (N^ - iy^[N^K,{Di) - n] = n{N^ - s)[{N^ - 

= J. On the other hand, Xu [22] showed that Ji{Di) > J, with 
equality if and only if Di is an OA{Ni,n, s,l). Therefore, Di must be an 
OA{Ni,n,s,l). 

(ii) Since the number of coincidences, 6ij{Di), must be an integer, it is 
easy to verify that, given Ji(Di) = J, J2{Di) achieves the minimum value 
of + ry(l — r/) when all 6ij{Di), i < j, take on only one of the two values, 
[JJ and [JJ + 1, where [x\ is the largest integer that does not exceed x. 
Then the result follows from (A.2). 

(iii) By Lemma 3, i^2,i(£>) = Lp_iK2,o(^) + 3-^X2(01). The result fol- 
lows from (ii) and the fact that K2fl{D) = s~'^n{n + s — 1). □ 

Proof of Lemma 5. By the binomial theorem, 

N N 

i=lj=l 

N N t . . 
i=lj=lr=0 ^ ^ 

By (8), (15) and (16), 

t , N 

KtiD^n) = N~'Yl I EEfe(^T)]*-^-L^-i 

+ ^"'E(0E E T.i^^ADT)Y-^{L;-L;^,) 

r=0 ^ ^ m=li£B,^jeB^ 
r=0 ^ ^ r=0 ^ ^ m=l 

= Ki,o(I?) + E ( I ) ^^i^'*-.o(^) 
+ ^"'E(I)(^p-^p-i)^*-'-(^i)- 

r=l ^ ' 
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Then the result fohows from Lemma 3 with some algebra. □ 

Proof of Lemma 6. First Asfi{D), A2,i{D) and A^^olD) are mini- 
mized sequentially because D has MA with respect to Wscf- Next, among 
designs with the same values of ^3^o(^)) A2,i{D) and A4^q{D), A^^q^D) is 
minimized because D has MA with respect to Wi, and A3^i{D) is minimized 
because D has MA with respect to Wscf - Continuing this type of argument 
shows that D has MA with respect to W2- □ 

Proof of Lemma 7. First 3^3,o(^) + ^2,i(^) = (3 - a)A3fl{D) + 
[aA3fi{D) + ^2,1(1?)] is minimized because both A3fi{D) and aA3fi{D) + 
^2,i{D) are minimized. For designs with the same value of 3A3^q{D) + 
^2,i{D), they must have the same values of ^3^o(^) and A2^i{D). Then 
^4,0 is minimized among designs with the minimum of 3A3fi{D) + A2^i{D) 
because D has MA with respect to Wgd- Among designs with the same values 
of 3A3^q{D) + A2^i{D) and A4^q{D), A^^q{D) is minimized because D has MA 
with respect to W2 and ^3^i(Z)) is minimized because D has MA with respect 
to Wsci] therefore, WA^^qID) + A^^ilD) is also minimized. For designs with 
the same values of 3^3,o(^) + A2,'i(-D), A^^iD) and 10^5,o(-D) +^3,1(1?), 
they must have the same A^fi{D) and ^3^1(1?) values. Continuing this type 
of argument shows that D has MA with respect to Wcc- D 

Proof of Theorem 3. Note that D achieves the lower bound in The- 
orem 2 if and only if Di achieves the lower bound in Lemma 4(ii). When the 
latter lower bound is achieved, K2{Di) is minimized and Kt{Di) is uniquely 
determined for t>3. By Lemma 3, Kt^i{D) is determined by Ktfl{D) for 
t > 3. Because Dt has MA, by Lemma 2, i^3,o(-D), i^^4,o(-D), . . . ,Knfi{D) are 
minimized sequentially. Then any combined sequence of (K3^o{D), K4fl(D), 
Knfi{D)) and (i^2,i(-D), K-i^i{D), Kn,i{D)) is also minimized sequen- 
tially as long as Kt^i{D) is minimized after Ktfi{D) for t = 2, . . . ,n. Hence, 
D has minimum moment aberration and MA with respect to Wsd, Wi and 
W2- Finally, because A2^i{D) is minimized among all possible designs, by 
Lemma 7, D has MA with respect to Wcc- □ 

Proof of Theorem 4. Necessity. When p = n — k — l,hy Lemma 4(i), 
the principal block Di is an OA{s, n, s, 1). Then it must contain a row of all 
zeros and other s — 1 rows without zeros. 

Sufficiency. Let u = {ui, . . . ,Un) be a row vector of Dt without zeros. 
Because none of Ui is zero, the linear equation Yl^Zi ^i^i = has s^~^~^ 
solutions over GF{s). Let F be an {n — k) x L^-k-i matrix, where the 
columns correspond to the solutions with the first nonzero element being 
unity. Clearly, F has rank n — k — 1, and it is an [n — k — 2)-flat in PG{n — 
k — l,s). On the other hand, Dt is an [n, n — k] linear code over GF{s). Let 
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T = (tij) be the {n — k) x n generator matrix of Dt- We need to show that 
T and F have no columns in common so that the resulting blocked design 
D = {Dt, Dp) is an RME design. 

Without loss of generality, let T = [In-k,E], where In-k is the n — k 
identity matrix and E is an {n — k)xk matrix. Because the row vectors of T 
form a basis for Dt, u can be uniquely represented as a linear combination of 
the row vectors of T. Then it is clear that ^27=1 ^ijUi = 7^ for j = 1, . . . , n. 
This proves that T and F have no columns in common. □ 

Proof of Theorem 5. We only need to prove the MA optimality. Fol- 
lowing the argument preceding the theorem, we can write an (s"'~^ : s"~'^~^) 
RME design as D = {Dt,Df), where Dt is a projection design of H^-k- 
Recall that the principal block Di is an OA{s,n,s,l); therefore, each level 
appears exactly once in each column. It is evident that 5ij{Di) = when i 7^ j 
and 6ij{Di) = n when i = j; hence, Kt{Di) = s~^'n!' for t > 0. By Lemma 3, 
Kt,i{D) is determined by Ktfl{D) for t > 0. By (11) and (12), At,i{D) is de- 
termined by Aifi{D), . . . , Atfi{D) uniquely. Thus, to sequentially minimize 
the sequences in (1), (2), (3) and (4), it is sufficient to sequentially minimize 
Aifi{D), . . . ,Anfl{D). Then the result follows from the condition that Dt 
has MA among all projection designs of Hn-k- D 

Proof of Theorem 6. Given K^fl{D), Kifl{D), . . . ,Kt,o{D), 3 <t < 
n, by Lemma 5, sequentially minimizing K2^i{D),K3^i{D), . . . , Kt-i^i{D) is 
equivalent to sequentially minimizing K3{D-an), K4^{D^^ri), . . . , Kt{D^^). Be- 
cause Dt has MA, by Lemma 2, Ks^Q{D),K/^fi(D), . . . ,Knfl{D) are mini- 
mized sequentially. Because D^^ has MA, K3{D^^^), Ki{D^ri), ■ ■ ■ , Kn{D^^^) 
are minimized sequentially. Then any combined sequence of {Ksfi{D) , Kifi{D) , 
. . . ,Knfi{D)) and {K2^i{D), K3^i{D), . . . ,Kn,i{D)) is also minimized sequen- 
tially as long as Kt-i^i{D) is minimized after Ktfl{D) for i = 3, . . . , n. Hence, 
D has minimum moment aberration and MA with respect to T^scf) and 
W2- By Lemma 5, K^{Dyia) = Kzfi^D) + 3K2,i(D) + constant; therefore, 
i^3,o(^) + 3^2,1 (^) is minimized and, by (13) and (14), A^fl{D) + A2,i{D) 
is minimized. Then, by Lemma 7, D has MA with respect to Wcc- CH 

Proof of Theorem 7. Because Dt has MA, A:i^Q{D) , Aifi{D) , . . . , 
Anfl{D) are minimized sequentially. Note that A2^i{D), A'i^i{D), and so on 
are minimized in Wi or W2 no sooner than in Wscf • Therefore, if D has MA 
with respect to Wsd, it must have MA with respect to Wi and W2. The MA 
Wcc optimality follows from Lemma 7. □ 
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